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Abstract: Different from the system in classical mathematics, a Smarandache system is 
a contradictory system in which an axiom behaves in at least two different ways within the 
same system, i.e., validated and invalided, or only invalided but in multiple distinct ways. 
Such systems exist extensively in the world, particularly, in our daily life. In this paper, we 
discuss such a kind of Smarandache system, i.e., non-solvable ordinary differential equation 
systems by a combinatorial approach, classify these systems and characterize their behaviors, 
particularly, the global stability, such as those of sum-stability and prod-stability of such 
linear and non-linear differential equations. Some applications of such systems to other 
sciences, such as those of globally controlling of infectious diseases, establishing dynamical 
equations of instable structure, particularly, the n-body problem and understanding global 


stability of matters with multilateral properties can be also found. 
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§1. Introduction 


Finding the exact solution of an equation system is a main but a difficult objective unless some 
special cases in classical mathematics. Contrary to this fact, what is about the non-solvable 
case for an equation system? In fact, such an equation system is nothing but a contradictory 
system, and characterized only by having no solution as a conclusion. But our world is overlap 
and hybrid. The number of non-solvable equations is much more than that of the solvable 
and such equation systems can be also applied for characterizing the behavior of things, which 
reflect the real appearances of things by that their complexity in our world. It should be noted 
that such non-solvable linear algebraic equation systems have been characterized recently by 
the author in the reference [7]. The main purpose of this paper is to characterize the behavior 


of such non-solvable ordinary differential equation systems. 
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Assume m, n > 1 to be integers in this paper. Let 
X = F(X) (DES") 
be an autonomous differential equation with F :R” — R” and F(0) = 0, particularly, let 
X=AX (LDES') 
be a linear differential equation system and 
a™) + aa +...4a,¢ =0 (LDE”) 


a linear differential equation of order n with 


G11 G12 +: Ain r1(t) filt, X) 
stats t t,X 

AeaM er ie S| OE | ass F(t, X) = falt,X) 
an aAn2 ea ann Ln(t) frit, X) 


where all a;, aij, 1 < 2,7 <n are real numbers with 
. . : - \T 
X = (41, 42,-++ ,¥n) 


and f;(t) is a continuous function on an interval [a,b] for integers 0 < i < n. The following 
result is well-known for the solutions of (LDES') and (LDE”) in references. 


Theorem 1.1([13]) Jf F(X) is continuous in 
U(X): |t-to] <a, |X — Xol] <b (@>0, b> 0) 


then there exists a solution X(t) of differential equation (DES') in the interval |t — to| < h, 
where h = min{a,b/M}, M = max || F(t, X)||- 


(t,X)€U (to, Xo) 
Theorem 1.2([13]) Let A; be the k;-fold zero of the characteristic equation 
det(A — AInxn) = |A- AInxn| = 0 
or the characteristic equation 


A” + aA" ++ +an_-1A + an =0 


with ky +ko+-+-:+k, =n. Then the general solution of (LDES') is 


n 
Se cB (the, 
i=1 
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where, c; is a constant, 3;,(t) is an n-dimensional vector consisting of polynomials in t deter- 
mined as follows 


ti 
3B ta1 
By (t) = 
tni 
tyit ar ty2 
3 tot +t 
Bo(t) = 21 22 
tnit + tn2 
re Tee Ae aoe Pe iegeg ral tik, 
kiy-1 t ky 2 
Z(t) = | ene + a ht tak 
it 
Tet | + gt? ae 
| t1(k1 41) ] 
3B t 
Brailt) = ae 
tn (ki +1) 
tyit ah ty2 
a tort +t 
Bri ro(t) = ae 
tnit + tn2 
Horrkst) phs—l 4 Um=kst2) phs—2 4... 4 ¢ 
(ks—1)! c= 7 
B Sah the ve winnie?) tks —2 tan 
By (t) = (ks—D) Gs 


t t 
n(n—ks+1) pky—1 n(n—kst2)pk,-2 1 
(k.—1)! t + a (ee) t + + thn 


with each ti; a real number for 1 <i,7 <n such that det((tij],,.,) A 0; 


Mt, tf l<i<kh; 
A2, tf ky +1 <i< ky; 


ay = 


Xs; if ky ko t---tkhs-p tl <ic<n. 


4 Linfan Mao 


The general solution of linear differential equation (LDE”) is 


s 


» (aa + cia? siete Ci(k;—1)t + Cik, er", 
i=1 


with constants qj, 1<i<s,l<j<ck. 


Such a vector family G,(t)e%*, 1 <i <n of the differential equation system (L DES") and 
a family t'e*!, 1 <1 < kj,1 <i <-s of the linear differential equation (LDE") are called the 
solution basis, denoted by 


B={ Bie |1<i<n} o CH{te*|1<i<s,1<1< k; }. 


We only consider autonomous differential systems in this paper. Theorem 1.2 implies that 
any linear differential equation system (LDES') of first order and any differential equation 
(LDE”) of order n with real coefficients are solvable. Thus a linear differential equation system 
of first order is non-solvable only if the number of equations is more than that of variables, and 
a differential equation system of order n > 2 is non-solvable only if the number of equations 
is more than 2. Generally, such a contradictory system, i.e., a Smarandache system [4]-[6] is 
defined following. 


Definition 1.3((4]-[6]) A rule R in a mathematical system (U;R) is said to be Smarandachely 
denied if it behaves in at least two different ways within the same set %, 1.e., validated and 
invalided, or only invalided but in multiple distinct ways. 

A Smarandache system (U;R) is a mathematical system which has at least one Smaran- 
dachely denied rule R. 


Generally, let (51;R1) (H2;Re2), +--+ ,(2mjRm) be mathematical systems, where R; is a 
rule on 4; for integers 1 <i < m. If for two integers i,j, 1 < i,j <m, 4; A 4; or Ly = 4; but 


Ri # Rj, then they are said to be different, otherwise, identical. We also aoe the conception 
of Smarandache multi-space defined following. 


Definition 1.4([4]-[6]) Let (X%1;Ri), (H2; Re), ++, (Um; Rm) be m > 2 mathematical spaces, 


different two by two. A Smarandache multi-space ¥ is a union U y; with rules R = U Ri on 
w=1 


= i.e., the rule Ry on %; for integers 1<1<m, denoted by (5: f). 


A Smarandache multi-space (5; R) inherits a combinatorial structure, i.e., a vertex-edge 


labeled graph defined following. 


Definition 1.5((4]-[6]) Let (E:R) be a Smarandache multi-space with S = US; and R = 
i=1 
U R;. Its underlying graph G [5,8 is a labeled simple graph defined by 


i=l 


V (G [=, R]) = {21, Ue,*--, Em}, 


B(G|5,&]) = { (4,2,) | LN D; 49,1< 4,5 <m} 
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with an edge labeling 


|F : (dy, U5) EE (c [S. R|) => ha O37 30) =W (B25) 5 
where w is a characteristic on 4; (|X; such that U,()X; is isomorphic to Xp) Xi if and only 
if @(Xi (U5) = w (Vel) X1) for integers 1 < i,j, k,l <m. 


Now for integers m, n > 1, let 


be a differential equation system with continuous F; : R” > R” such that F;(0) = 0, particu- 
larly, let 
X= A,X,-++,X =Az,X,---,X = AnX (LDES},) 


be a linear ordinary differential equation system of first order and 


gh) 4 al) g(m—1) vee ale =0 
a) +ablem-) 4...4 022 =0 
: . (LDE®) 


g®™ ol Es all gm) + ners oe alla = 0 


a linear differential equation system of order n with 


k k k 
Ge ale aeceale #1(t) 
k k k 

ice: ay a On | gy | tale) 
re ee alt) an(t) 


[A] 


where each a;;° is a real number for integers 0 << k <m, 1<i,j <n. 


Definition 1.6 An ordinary differential equation system (DES},) or (LDES},) (or (LDE™)) 
are called non-solvable if there are no function X(t) (or x(t)) hold with (DES},) or (LDES),) 
(or (LDE”.)) unless the constants. 


The main purpose of this paper is to find contradictory ordinary differential equation 
systems, characterize the non-solvable spaces of such differential equation systems. For such 
objective, we are needed to extend the conception of solution of linear differential equations in 


classical mathematics following. 


Definition 1.7 Let S° be the solution basis of the ith equation in (DES},). The V-solvable, A- 
solvable and non-solvable spaces of differential equation system (DES) are respectively defined 
by 


Ls? (5? and (Js? 1) 5°, 
t=1 t=1 t=1 t=1 


where S° is the solution space of the ith equation in (DES}). 
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According to Theorem 1.2, the general solution of the ith differential equation in (LDES} ) 
or the ith differential equation system in (ZLDE?) is a linear space spanned by the elements 
in the solution basis A; or G; for integers 1 < i < m. Thus we can simplify the vertex-edge 
labeled graph G [> Rl replaced each }7; by the solution basis 4; (or @) and >7;()>0, by 
BiB; (or 6 ()G;) if A VB; AO (or GG; AO) for integers 1 < i, 7 < m. Such a vertex- 
edge labeled graph is called the basis graph of (LDES},) ((LDE”,)), denoted respectively by 
G[LDES}] or G[LDE®] and the underlying graph of G[LDES}] or G[LDE™], i.e., cleared 
away all labels on G[LDES1] or G[LDE®,] are denoted by G[LDES}] or G[LDE®). 

Notice that a oo = U S°, i.e., the non-solvable space is empty only if m = 1 in 


w=1 i=l 
(LDEq). Thus G[LDES'| ~ Kk, or G[LDE"| ~ kK, only if m = 1. But in general, the 


basis graph G[LDES},] or G[LDE®™] is not trivial. For example, let m = 4 and 4? = 
{errt, er2t ersth £3 = {erst erat erst}. BS = {errt erst erst} and BY = {erst erst ret}. 
where \;, 1 < i < 6 are real numbers different two by two. Then its edge-labeled graph 
G[LDES}] or G[LDE™] is shown in Fig.1.1. 


as 
@) Be’ 


{ert erst} {erat ersth 


(Ze 


Fig.1.1 


If some functions F;(X), 1 < i < m are non-linear in (DES},), we can linearize these 
non-linear equations X = F;(X) at the point 0, ie., if 


F(X) = F/(0)X + Ri(X), 


where F/(0) is an n x n matrix, we replace the ith equation X = F;(X) by a linear differential 


equation 


X= F/O)X 


in (DES},). Whence, we get a uniquely linear differential equation system (LDES},) from 
(DES}.) and its basis graph G[LDES!]. Such a basis graph G[LDES}] of linearized differen- 
tial equation system (DES},) is defined to be the linearized basis graph of (DES},) and denoted 
by G[DES}]. 

All of these notions will contribute to the characterizing of non-solvable differential equation 
systems. For terminologies and notations not mentioned here, we follow the [13] for differential 
equations, [2] for linear algebra, [3]-[6], [11]-[12] for graphs and Smarandache systems, and [1], 
[12] for mechanics. 
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§2. Non-Solvable Linear Ordinary Differential Equations 


2.1 Characteristics of Non-Solvable Linear Ordinary Differential Equations 


First, we know the following conclusion for non-solvable linear differential equation systems 
(LDES}) or (LDE®). 


Theorem 2.1 The differential equation system (LDES},) is solvable if and only if 
(|A1 — AInxns |A2 — Mnxn|s+++ 5 |Am — Inxn|) #1 
i.e., (LDEgq) is non-solvable if and only if 
(|A1 — AInxn,|A2 — Inxnl,+ ++ 5 [Am — AInxn|) = 1. 
Similarly, the differential equation system (LDE®) is solvable if and only if 
(Pi(A), PaQa),-++ s Pm(A)) 4 1, 
i.e., (LDE?) is non-solvable if and only if 


(Pid), Pa(A), +++ > Pm(A)) = 1, 


0] 


where P;(X) = A” + alt] n-1 very El » for integers 1<i<m. 


[ 
i(n—1)> +4; 


Proof Let Ai1, Ai2,-++ , Ain be the n solutions of equation |A; — AInxn| = 0 and &; the 
solution basis of ith differential equation in (LDES},) or (LDE”,) for integers 1 < i < m. 
Clearly, if (LDES!,) ((LDE®.)) is solvable, then 


() Bi #9, Le., (Prins diz. rin AO 
t=1 t=1 


by Definition 1.5 and Theorem 1.2. Choose Ao € ial {\i1,Ai25°°+ Ain}. Then (A — Xo) is a 


common divisor of these polynomials |A1 — AInxn, oo Mnxn|)+++ )|Am — AInxn|. Thus 
(\Ai — AInxn; |A2 — AInxnls-++ y|Am — AInxn|) 4 1. 
Conversely, if 
(|Ai — AInxn,|A2 — AInxn|y+++ s|Am — Alnxnl) 4 1, 


let (A—Ao1); (A—Ao2),°++ » (A—Aor) be all the common divisors of polynomials |A,—AInxn, |A2—- 
Mnxnl,-++ ,|Am — MInxn|, where Ao; A Aoz if i A 7 for 1 <i, 7 <l. Then it is clear that 


Cye*? + Coe*? +--+. + Chero 


is a solution of (LEDq) ((LDE”.)) for constants Cy,C2,--- ,C). 


For discussing the non-solvable space of a linear differential equation system (LEDS1,) or 
(LDE”.) in details, we introduce the following conception. 
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Definition 2.2 For two integers 1< i,j <m, the differential equations 


dx; 
AX 
dt 
1 
(LDES}.) 
—i=A,X 
dt q 
in (LDES}) or 
20 + als) 4. pale =o i 
B00 all gn) ss alle _~90 ” 


in (LDE”.) are parallel if Bi(.\B; =9. 
Then, the following conclusion is clear. 


Theorem 2.3 For two integers 1 < i,j <_m, two differential equations (LDES};) (or (LDE}.)) 
are parallel if and only if 


(|Ai] — nxn; |Agl — nxn) = 1 (or (Pi(A), PQA) = 1), 


where (f(x), g(x)) is the least common divisor of f(x) and g(x), Py(A) =A" + al’) AME eee 


ath ay’ + all fork =1,9. 


Proof By definition, two differential equations (LEDSj,) in (LDESj,) are parallel if and 
only if the characteristic equations 


|A; — Mnxn|=0 and |A; —Anxn| =0 


have no same roots. Thus the polynomials |A;| — AInxn and |A;| — AInxn are coprime, which 
means that 


(\Ai — Alnxn, |Aj — Alnxn) = 1. 


Similarly, two differential equations (LED7,) in (LDE7,) are parallel if and only if the 
characteristic equations P;(A) = 0 and P;(A) = 0 have no same roots, i.e., (P;(A), Pj(A)) = 1. 


Let f(x) =agz™ +4,27 1 +--+ +Gm12+4m, g(x) = box” + bya" +--+ + bya t by 
with roots 71,%2,-++,%m and y1, y2,°-: , Yn, respectively. A resultant R(f,g) of f(x) and g(x) 
is defined by 


R(f,9) = 45°05 [ [ (ei — ys). 
tJ 
The following result is well-known in polynomial algebra. 


Theorem 2.4 Let f(x) = agxz™ +a,2™ 1 +--+ +@m-1t+4m, g(x) = box” + bya" 1 Hee + 
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by—-12 + by with roots 71,%2,°°: ,Lm and y1,Y2,°°* ,Yn, respectively. Define a matrix 
ao ay eee Am 0 ae 0 0 
0 ao ay Bie ie Am 0 eee 0 
0 0 0 ao ay am 
V(f,g) = 
bile os Be GY sae 
Or Be Be a2 Dy ase, NO 
0 0 0 bo by Dn 
Then 


R(f, 9) = detV(f, g)- 


We get the following result immediately by Theorem 2.3. 


Corollary 2.5 (1) For two integers 1 < i,j < m, two differential equations (LDES;;) are 
parallel in (LDES},) if and only if 


R(JAi — Mnxnl,|Aj — nxn) #9; 
particularly, the homogenous equations 
V(\Ai — Inxn|, Aj — Lnxnl)X = 0 
have only solution (0,0,---,0)" if |A; — AInxn| = @oA” + aA"! + ++» + Gn-1A + Gn and 
a 


|A; — AInxn| = boA” + bart +++ + bn-iA + bn. 
(2) For two integers 1 < i,j < m, two differential equations (LDE};) are parallel in 
(LDE”) if and only if 
R(Pi(A), Py(A)) 0, 
particularly, the homogenous equations V(P;(A), Pj(A))X =0 have only solution (0,0,-++ ,0)7. 
ola 


Proof Clearly, |Ai — AInxn| and |A; — AInxn| have no same roots if and only if 
R(|Ai — Mnxn|,|Aj — Inxn|) #9, 


which implies that the two differential equations (LEDS/,) are parallel in (LEDS;,,) and the 
homogenous equations 
VAT eal (xe 6 


have only solution (0,0,--- ,0)7. That is the conclusion (1). The proof for the conclusion (2) 
— 


2n 


is similar. 


10 Linfan Mao 


Applying Corollary 2.5, we can determine that an edge (4, Z;) does not exist in G[LDES, | 
or G[LDE*,] if and only if the ith differential equation is parallel with the jth differential equa- 
tion in (LDES},) or (LDE”,). This fact enables one to know the following result on linear 


non-solvable differential equation systems. 


Corollary 2.6 A linear differential equation system (LDES},) or (LDE®,) is non-solvable if 
G(LDES!,) % Km or G(LDE®,) & Km for integers m,n > 1. 


2.2. A Combinatorial Classification of Linear Differential Equations 


There is a natural relation between linear differential equations and basis graphs shown in the 


following result. 


Theorem 2.7 Every linear homogeneous differential equation system (LDES),) (or (LDE™)) 
uniquely determines a basis graph GILDES},] (G[LDE™]) inherited in (LDES}) (or in(LDE”)). 
Conversely, every basis graph G uniquely determines a homogeneous differential equation system 


(LDES}) ( or (LDE”)) such that G[LDES},] ~ G (or G[LDE”] ~ G). 


Proof By Definition 1.4, every linear homogeneous differential equation system (LDES} ) 
or (LDE”.) inherits a basis graph G[LDES}] or G[LDE®], which is uniquely determined by 
(LDES},) or (LDE”). 

Now let G be a basis graph. For Vu € V(G), let the basis Z, at the vertex v be A, = 
{ B,(t)ev* | 1 <i< ny} with 


A, tif 1St<sky; 


A2, tf ky +1 <i< ko; 
a= 


We construct a linear homogeneous differential equation (LDES') associated at the vertex v. 
By Theorem 1.2, we know the matrix 


tir ota ts tn, 
T= to1 toa +++ tan, 
bad tn,2 rae? Cie 


is non-degenerate. For an integer i, 1<1i< s, let 
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be a Jordan black of k; x k; and 


— = AX (LDES') 


with X = [21 (t), v2(t),--- ,an,(t)|” is nothing but 4, by Theorem 1.2. Notice that the Jordan 
black and the matrix T are uniquely determined by &,. Thus the linear homogeneous differen- 
tial equation (L DES") is uniquely determined by &,. It should be noted that this construction 
can be processed on each vertex v € V(G). We finally get a linear homogeneous differential 
equation system (LDES1},), which is uniquely determined by the basis graph G. 

Similarly, we construct the linear homogeneous differential equation system (ZLDE”) for 
the basis graph G. In fact, for Vu € V(G), let the basis Z, at the vertex u be B, = { tle®* | 1 < 
i<s,1<Il<k;}. Notice that A; should be a k;-fold zero of the characteristic equation P(A) = 0 
with ky +ko +--+: +k, =n. Thus P(Ay) = P’(Ay) =-+: = PM -D(A;) = 0 but P@)(A;) 4 0 
for integers 1 <i < s. Define a polynomial P,,(A) following 


s 


P,Q) = [JQ—>)" 
i=l 
associated with the vertex u. Let its expansion be 
Pu(A) =A" + gx a a Qu(n—1)A + Gun- 
Now we construct a linear homogeneous differential equation 


gh”) + anja) peek Ht Qu(n—1)e" + Adynt = 0 (LP DES) 


associated with the vertex u. Then by Theorem 1.2 we know that the basis solution of (L DE”) 
is just G,. Notices that such a linear homogeneous differential equation (LDE”) is uniquely 


constructed. Processing this construction for every vertex u € V(G), we get a linear homoge- 


neous differential equation system (LDE”,). This completes the proof. 
Example 2.8 Let (LDE*) be the following linear homogeneous differential equation system 


&— 34+ 2x = 0 1 
z—5x%+ 6x = 0 2 
&—7z+12¢% =0 3 
X— 9x + 20x = 0 4 
z—11la+ 307% = 0 5 
xz—7x+6x = 0 6 


Sst Ne US et RE en 
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2 


d d 
where % = <* and ¢ = —. Then the solution basis of equations (1) — (6) are respectively 


fete}, {e7# eft, fer eM}, fet, eF4}, {e54, e%1, {e% e4} and its basis graph is shown in 
Fig.2.1. 
{ete} fe} {e2t, ety 
{e"} {e**} 
{et et} {et At} 
{et} {eAry 
{et eSt} [e51} fet et} 


Fig.2.1 The basis graph H 


Theorem 2.7 enables one to extend the conception of solution of linear differential equation 
to the following. 


Definition 2.9 A basis graph G[LDES},] (or G[LDE®]) is called the graph solution of the 
linear homogeneous differential equation system (LDES},) (or (LDE™)), abbreviated to G- 


solution. 


The following result is an immediately conclusion of Theorem 3.1 by definition. 


Theorem 2.10 Every linear homogeneous differential equation system (LDES},) (or (LDE™)) 
has a unique G-solution, and for every basis graph H, there is a unique linear homogeneous 
differential equation system (LDES},) (or (LDE®,)) with G-solution H. 


Theorem 2.10 implies that one can classifies the linear homogeneous differential equation 
systems by those of basis graphs. 


Definition 2.11 Let (LDES},), (LDES},)’ (or (LDE™), (LDE”)') be two linear homo- 
geneous differential equation systems with G-solutions H, H’. They are called combinato- 
rially equivalent if there is an isomorphism yp : H — H’, thus there is an isomorphism 
yp: H — H' of graph and labelings 0, 7 on H and H' respectively such that p(x) = Ty(x) for 
Vc € V(H)J E(H), denoted by (LDES},) & (LDES1)! (or (LDE”) © (LDE®)’). 


{e he any {e atte 2t e abt 
e 34} 


fe- st, ee 


goat 


{e Ste 6t\ fe Bel de Ate ae 


Fig.2.2. The basis graph H’ 
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Example 2.12 Let (LDE?.)' be the following linear homogeneous differential equation system 


+34 +2¢ =0 1) 
#+5é+6r =0 2) 
#+7¢+12c =0 3) 
#494 +202 =0 4) 
#+11¢+302=0 (5) 
&+74 + 6c =0 6) 


Then its basis graph is shown in Fig.2.2. 
Let y: H — H’ be determined by y({e>*, e%7*}) = {e-*#, es") and 


v({er*, ert} (emery = fe ie fe heh 


for integers 1 < i,k <6 and j =i+1 = 6(mod6), 1 = k+1= 6(mod6). Then it is clear that 
H & H'. Thus (LDE”)' is combinatorially equivalent to the linear homogeneous differential 
equation system (LDE?.) appeared in Example 2.8. 


Definition 2.13 Let G be a simple graph. A vertex-edge labeled graph 6: G— Z* is called 
integral if 0(uv) < min{@(u), O(v)} for Vuv € E(G), denoted by G. 

Let Ge and Gy be two integral labeled graphs. They are called identical if G1 x G2 and 
O(x) = r(y(x)) for any graph isomorphism y and Vx € V(G1) UE(Gi), denoted by Gi? = GP. 


For example, these labeled graphs shown in Fig.2.3 are all integral on K4—e, but Ot, = ene 
GREG? 
1 # G3. 


3 2 4 4 a 3 3 1 3 
1 2 2 1 2 2 
4 2 3 3 4 4 1 4 
Ce Gr Gy 
Fig.2.3 


Let G[LDES}] (G[LDE”]) be a basis graph of the linear homogeneous differential equa- 
tion system (LDES},) (or (LDE™)) labeled each v € V(G[LDES}]) (or v € V(G[LDE™])) 
by &,. We are easily get a vertex-edge labeled graph by relabeling v € V(G[LDES}]) (or 
v € V(G[LDE™])) by |#| and wv € E(G[LDES}]) (or uv € E(G[LDE™])) by |Buf\ Fol. 
Obviously, such a vertex-edge labeled graph is integral, and denoted by G/[LDES},] (or 
G'|LDE®™]). The following result completely characterizes combinatorially equivalent linear 
homogeneous differential equation systems. 


Theorem 2.14 Let (LDES},), (LDES},)' (or (LDE®,), (LDE”)') be two linear homogeneous 
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differential equation systems with integral labeled graphs H, H’. Then (LDES},) ag (LDES})! 


(or (LDE”) & (LDE®)') if and only if H = H’. 


Proof Clearly, H = H' if (LDES},) & (LDES},)' (or (LDE”.) © (LDE”.)’) by defini- 
tion. We prove the converse, i.e., if H = H’ then there must be (LDES?,) © (LDES1,)’ (or 
(LDE”) & (LDE®)’). 

Notice that there is an objection between two finite sets 5), So if and only if |S,| = |So|. 
Let tT be a1 —1 mapping from ¥, on basis graph G[LDES})] (or basis graph G[LDE™]) to 
By on basis graph G[LDES} |’ (or basis graph G[LDE™]') for v,v' € V(H’). Now if H = H’, 
we can easily extend the identical isomorphism zdy on graph H to a 1— 1 mapping id}, : 
G[LDES}] — G[LDES}]' (or id}, : G[LDE®™] — G[LDE™]’) with labelings 0: v — &, and 

i. v0 + By on G[LDES}], G[LDES}]' (or basis graphs G[LDE™], G[LDE™]'). Then 
it is an immediately to check that id%,0(x) = 6’r(x) for Vz € V(G[LDES) |) E(G[LDES,]) 
(or for Vz € V(G[LDE?,])  E(G[LDE”,])). Thus id}, is an isomorphism between basis graphs 
G[LDES}] and G[LDES}]' (or GILDE®] and G[LDE]'). Thus (LDES}) ‘“' (LDES},y 


id’ 
(or (LDE”) ~ (LDE”)’). This completes the proof. 


According to Theorem 2.14, all linear homogeneous differential equation systems (L DES} ) 
or (LDE?.) can be classified by G-solutions into the following classes: 


Class 1. G[LDES!] ~ Km or G[LDE®] ~ Ky for integers m,n > 1. 


The G-solutions of differential equation systems are labeled by solution bases on Kj, and 
any two linear differential equations in (LDES},) or (LDE”.) are parallel, which characterizes 


m isolated systems in this class. 


For example, the following differential equation system 


i+ 3x + 2x =0 
«—5x+ 6x = 0 
@+ 2¢ —- 3x =0 


is of Class 1. 
Class 2. G[LDES1] ~ Km or GILDE®] ~ Ky, for integers m,n > 1. 


The G-solutions of differential equation systems are labeled by solution bases on complete 
graphs K, in this class. By Corollary 2.6, we know that G[LDES!] ~ Km or GILDE"| ~ Km, 
if (LDES},) or (LDE”.) is solvable. In fact, this implies that 


() B= [) uf) a) 46. 


vEV (Km) u,vEV (Km) 


Otherwise, (LDES},) or (LDE”,) is non-solvable. 


For example, the underlying graphs of linear differential equation systems (A) and (B) in 
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the following 


Z-— 32+ 27 =0 
Z-— 34 +27 =0 
r—-x=0 


(A) (B) # — 5a +62 =0 
% —4¢ +32 =0 
# —4¢ +32 =0 


i+ 2x — 3x =0 


are respectively K4, K3. It is easily to know that (A) is solvable, but (B) is not. 


Class 3. G[LDES!] ~ G or G[LDE®] ~ G with |G| =m but G # Km, Km for integers 
m,n > 1. 


The G-solutions of differential equation systems are labeled by solution bases on G and all 
linear differential equation systems (LDES},) or (LDE®,) are non-solvable in this class, such 
as those shown in Example 2.12. 


2.3 Global Stability of Linear Differential Equations 


The following result on the initial problem of (LDES') and (LDE”) are well-known for differ- 
ential equations. 


Lemma 2.15([13]) For t € [0,00), there is a unique solution X(t) for the linear homogeneous 
differential equation system 


“ = AX (L"DES") 


with X(0) = Xo and a unique solution for 
oh) + aya?) +---+ta,v=0 (L? DE”) 
with 2(0) = 29, 2'(0) = 2,--- ,a("-Y(0) = ag 


Applying Lemma 2.15, we get easily a conclusion on the G-solution of (LDES},) with 
X,(0) = X® for Vu € V(G) or (LDE®) with x(0) = x, 2'(0) = 2),--- ,a@-Y(0) = a") by 
Theorem 2.10 following. 


Theorem 2.16 For t © [0,00), there is a unique G-solution for a linear homogeneous dif- 
ferential equation systems (LDES},) with initial value X,(0) or (LDE™) with initial values 
ay(0),2/,(0),--- , a6") (0) for Vu € V(G). 


For discussing the stability of linear homogeneous differential equations, we introduce the 
conceptions of zero G-solution and equilibrium point of that (LDES},) or (LDE”,) following. 


Definition 2.17 A G-solution of a linear differential equation system (LDES},) with initial 
value X,(0) or (LDEP.) with initial values x,(0), x/,(0),--- ,af"— (0) for Vu € V(G) is called 
a zero G-solution if each label B; of G is replaced by (0,--- ,0) (|\Bi| times) and Bi (|B; by 
(0,--- ,0) (Ai A;| times) for integers 1<i,7 <m. 
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Definition 2.18 Let dX/dt = A,X, oO) + aye") +--+ aynt = 0 be differential equations 
associated with verter v and H a spanning subgraph of G[LDES},| (or G[LDE®]). A point 
X* € R” is called a H-equilibrium point if AyX* = 0 in (LDES},) with initial value X,(0) 
or (X*)" + yi (X*)"1 + +++ + Gyn X* = 0 in (LDE”,) with initial values x,(0), x/,(0),---, 
a" (0) for Vu € V(A). 


We consider only two kind of stabilities on the zero G-solution of linear homogeneous 


differential equations in this section. One is the sum-stability. Another is the prod-stability. 


2.3.1 Sum-Stability 


Definition 2.19 Let H be a spanning subgraph of GILDES},] or G[LDE™] of the linear 
homogeneous differential equation systems (LDES},) with initial value X,(0) or (LDE”) with 
initial values vy(0),2/,(0),--+, af Yo). Then G[LDES}] or GILDE®] is called sum-stable 
or asymptotically sum-stable on H if for all solutions Y,(t), v € V(H) of the linear differential 
equations of (LDES},) or (LDE™) with |Y,(0)—Xv(0)| < 5, exists for allt >0,| SS Y,(t)- 
veV (HA) 
So Xy(t)| <¢, or furthermore, lim| S> Y(t)- SS Xy(t)| =0. 
veV(H) $0 veV(H) veV(H) 
Clearly, an asymptotic sum-stability implies the sum-stability of that G[LDES},] or G[LDE” ]. 
The next result shows the relation of sum-stability with that of classical stability. 


Theorem 2.20 For a G-solution G[LDES} ] of (LDES},) with initial value X,,(0) (or GILDE™] 
of (LDEP.) with initial values x,(0),x/,(0),---, a" (0)), let H be a spanning subgraph of 
G[LDES}] (or G[LDE™]) and X* an equilibrium point on subgraphs H. If G[LDES}] (or 
G[LDE™]) is stable on any Vu € V(H), then GILDES}] (or G[LDE®]) is sum-stable on H. 
Furthermore, if G[LDES},] (or G[LDE®]) is asymptotically sum-stable for at least one vertex 
v€V(H), then GILDES}] (or G[LDE™]) is asymptotically sum-stable on H. 


Proof Notice that 
| SS pKO- SS pX Hl < So plv(t) - Xl) 
a (A) veEV(H) veEV(H) 


and 


lim| SY) pY(t)- DY) peXOlS DS) po lim|¥.() - Xv(@)]. 


veV (HA) veV(H) veV(H) 


Then the conclusion on sum-stability follows. 
For linear homogenous differential equations (LDES") (or (LDE")), the following result 
on stability of its solution X(t) = 0 (or x(t) = 0) is well-known. 


Lemma 2.21 Let y = max{ Red] |A— Alnxn| = 0}. Then the stability of the trivial solution 
X(t) = 0 of linear homogenous differential equations (LDES') (or x(t) = 0 of (LDE™)) is 


determined as follows: 


(1) ify <0, then it is asymptotically stable; 
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(2) ify > 0, then it is unstable; 
(3) of y = 0, then it is not asymptotically stable, and stable if and only if m'(A) = m(A) 
for every with ReX = 0, where m(A) is the algebraic multiplicity and m’(A) the dimension of 


eigenspace of X. 


By Theorem 2.20 and Lemma 2.21, the following result on the stability of zero G-solution 
of (LDES},) and (LDE”) is obtained. 


Theorem 2.22 A zero G-solution of linear homogenous differential equation systems (LDES} ) 
(or (LDE®.)) is asymptotically sum-stable on a spanning subgraph H of GILDES},] (or G[LDE”]) 
if and only if Rea, < 0 for each B,,(t)e%* € B, in (LDES') or Redy <0 for each t'-e* € ©, 
in (LDE?.) hold for Vu € V(A). 


Proof The sufficiency is an immediately conclusion of Theorem 2.20. 
Conversely, if there is a vertex v € V(H) such that Rea, > 0 for B,(t)e%* € Ay in 
(LDES') or Red, > 0 for t've*! € G, in (LDE®,), then we are easily knowing that 
jim B,(t)e** — co 


if a, > 0 or B,(t) Aconstant, and 


lim t' eet 5 00 
t—00 


if Ay > 0 or ly, > 0, which implies that the zero G-solution of linear homogenous differential 


equation systems (LDES") or (LDE") is not asymptotically sum-stable on H. 


The following result of Hurwitz on real number of eigenvalue of a characteristic polynomial 
is useful for determining the asymptotically stability of the zero G-solution of (LDES!},) and 
(LDE”). 


Lemma 2.23 Let P(X) =A" +a,A\"-1 +++» +an-1A + an be a polynomial with real coefficients 
aj, 1<i<nand 


at 1 0 0 
a3 ag ay, 0 0 
ay, 1 
Ai = lai|, Ag = ; An = a5 a4 a3 @2 ay 0 0 
a3 «a2 
0 Reet: Qn 


Then ReX < 0 for all roots X of P(A) if and only if A; > 0 for integers 1<i<n. 


Thus, we get the following result by Theorem 2.22 and lemma 2.23. 


Corollary 2.24 Let Aj, A$,--- ,A? be the associated determinants with characteristic polyno- 
mials determined in Lemma 4.8 for Vv € V(G[LDES}]) or V(G[LDE™]). Then for a spanning 
subgraph H < G[LDES!] or G[LDE™], the zero G-solutions of (LDES},) and (LDE™,) is 
asymptotically sum-stable on H if AY > 0,A3 >0,---,A® > 0 for Vu € V(A). 
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Particularly, if n = 2, we are easily knowing that Red < 0 for all roots A of P(A) if and 
only if a, > 0 and ag > 0 by Lemma 2.23. We get the following result. 


Corollary 2.25 Let H < G[LDES}] or G[LDE”] be a spanning subgraph. If the characteristic 
polynomials are \7 + a¥\+a3 for v € V(H) in (LDES},) (or (L"DE?,)), then the zero G- 
solutions of (LDES},) and (LDE?,) is asymptotically sum-stable on H if at > 0, a} > 0 for 
Vue V(A). 


2.3.2 Prod-Stability 


Definition 2.26 Let H be a spanning subgraph of G[LDES},] or G[LDE™] of the linear 
homogeneous differential equation systems (LDES},) with initial value X,(0) or (LDE®,) with 
initial values xy(0),2/,(0),---, a0). Then GILDES}] or G[LDE™] is called prod-stable 


or asymptotically prod-stable on H if for all solutions Y,(t), v € V(#A) of the linear differential 
equations of (LDES},) or (LDE™) with |Y,(0)—Xv(0)| < 5, exists for allt >0,| T] Y.(t)- 
veV (HA) 
II Xv(t)| <, or furthermore, lim| [] Y,(t)- J] Xv(t)| =0. 
veV(H) t>0  veV(H) veV(H) 
We know the following result on the prod-stability of linear differential equation system 
(LDES1) and (LDE"). 


Theorem 2.27 A zero G-solution of linear homogenous differential equation systems (LDES} ) 
(or (LDE™ )) is asymptotically prod-stable on a spanning subgraph H of GILDES},] (or G[LDE”]) 


if and only if S > Reay < 0 for each B,(t)e%t € By in (LDES') or SY > Redry <0 for 
veEV(H) veEV(H) 


each tert € &, in (LDE™). 


Proof Applying Theorem 1.2, we know that a solution X,(t) at the vertex v has the form 


Whence, 
II Xy (t) = ye cB, (t)er" 
veV(H) veV(H) t=1 
n _ n = XS at 
= | IL owe }=)> I] able 
i=1 veV(H) i=1 veV(H) 


Whence, the zero G-solution of homogenous (LDES?,) (or (LDE®.)) is asymptotically sum- 
stable on subgraph H if and only if >> Rea, < 0 for VG,(t)e%’ € BZ, in (LDES') or 
veEV (A) 
d> Red, <0 for Vtlvert € &, in (LDE™). 
veV (HA) 


Applying Theorem 2.22, the following conclusion is a corollary of Theorem 2.27. 


Corollary 2.28 A zero G-solution of linear homogenous differential equation systems (LDES}. ) 
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(or (LDE”)) is asymptotically prod-stable if it is asymptotically sum-stable on a spanning 


subgraph H of G[LDES}] (or G[LDE®]). Particularly, it is asymptotically prod-stable if the 


zero solution 0 is stable on Vu € V(A). 


Example 2.29 Let a G-solution of (LDES},) or (LDE”,) be the basis graph shown in Fig.2.4, 
where v; = {e724, e738}, uy = {e734 eM}, ug = fet 8, 84}, vy = fe75t e8t, eB}. 
vs = {e*,e ©}. ug = {e *,e 7, e *}. Then the zero G-solution is sum-stable on the triangle 


U4U5 U6, but it is not on the triangle v,v2v3. In fact, it is prod-stable on the triangle v,v2v3. 


ve  {e**} V1 


{e“} 
U5 tere | Aewt 


{e— 


V4 feel. Us 


Fig.2.4 A basis graph 


§3. Global Stability of Non-Solvable Non-Linear Differential Equations 


For differential equation system (DES},), we consider the stability of its zero G-solution of 
linearized differential equation system (LDES},) in this section. 


3.1 Global Stability of Non-Solvable Differential Equations 


Definition 3.1 Let H be a spanning subgraph of G[DES}] of the linearized differential equation 
systems (DES},) with initial value X,(0). A point X* € R” is called a H-equilibrium point of 
differential equation system (DES},) if f.(X*) =0 for Vu € V(A). 


Clearly, 0 is a H-equilibrium point for any spanning subgraph H of G[DES},] by definition. 
Whence, its zero G-solution of linearized differential equation system (LDES1}) is a solution 
of (DES}). 


Definition 3.2 Let H be a spanning subgraph of G[DES}] of the linearized differential equation 
systems (DES?) with initial value X,(0). Then G[DES}] is called sum-stable or asymptoti- 
cally sum-stable on H if for all solutions Y,(t), v € V(H) of (DES},) with ||Y,(0) — Xy(0)|| < dy 
exists for all t > 0, 


+ KwO- YL XO] <a 


veV(H) veV(H) 


or furthermore, 
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lim S> YwYO- So XM] =0, 


veV(H) veV(H) 
and prod-stable or asymptotically prod-stable on H if for all solutions Y,(t), v € V(H) of 
(DES1,) with ||Y,(0) — X,(0)|| < 6, exists for all t > 0, 


I] “®- TI xl <e 


veV(H) veV(H) 


or furthermore, 


im|| J] w@®- J] x] =o. 


t—0 
veV(H) veV(H) 


Clearly, the asymptotically sum-stability or prod-stability implies respectively that the 
sum-stability or prod-stability. 


Then we get the following result on the sum-stability and prod-stability of the zero G- 
solution of (DES}). 


Theorem 3.3 For a G-solution G[DES!},] of differential equation systems (DES},) with initial 
value X,,(0), let Hy, Hz be spanning subgraphs of G[DES}]. If the zero G-solution of (DES),) 
is sum-stable or asymptotically sum-stable on H, and Hz, then the zero G-solution of (DES},) 
is sum-stable or asymptotically sum-stable on Hy U Ho. 

Similarly, if the zero G-solution of (DES},) is prod-stable or asymptotically prod-stable on 
H, and X,(t) is bounded for Yv € V(Hz2), then the zero G-solution of (DES},) is prod-stable 
or asymptotically prod-stable on Hy |) Ho. 


Proof Notice that 
|X1 + X2|| < |[Xa|] + [|X|] and ||X1Xel] < || X1||||Xell 


in R”. We know that 


S- X,()|| = S> Xd ()+ SS Xe) 


veV (1) U V(A2) veV (1) veV (H2) 


IA 
M 
= 
+ 
M 
s 


and 


II X,(t) 


veEV(A1)U V(A2) veV (Ai) veEV(H2) 


IA ll 
Ss aS 
i. Le 
— 
| 
is 
x S 
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Whence, 


vE€V(A1)U V(A2) 


if € =€, +€2 with 


veEV(A1) 


or 


lim 
t—0 


veEV(A1) 


Se Xue) 
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<e€ or lim 
veEV(M1)U V(A2) 


<e& and 


S> Xp (t)]] < e2 


vEV (Ha) 


=0 and lim =0. 


> 


Ss) X(t) 


vEV (He) 


This is the conclusion (1). For the conclusion (2), notice that 


I 


vEV(A1)U V(A2) 


X,(t) 


I] xt 


veEV(H1) 


<} [[ x@l] IT] x.@]| < Me 
veV (1) veV (H2) 
)I| <e and I] x] <™. 


vEV (He) 


Consequently, the zero G-solution of (DES},) is prod-stable or asymptotically prod-stable on 


A, U Ap. 


Theorem 3.3 enables one to get 


the following conclusion which establishes the relation of 


stability of differential equations at vertices with that of sum-stability and prod-stability. 


Corollary 3.4 For a G-solution G[DES},] of differential equation system (DES},) with initial 
value X,(0), let H be a spanning subgraph of G[DES},]. If the zero solution is stable or 
asymptotically stable at each vertex v € V(A), then it is sum-stable, or asymptotically sum- 
stable and if the zero solution is stable or asymptotically stable in a vertex u € V(A) and X(t) 
is bounded for Vv € V(H) \ {u}, then it is prod-stable, or asymptotically prod-stable on H. 


It should be noted that the converse of Theorem 3.3 is not always true. For example, let 


S- X,(t)|| <a+e and S- Xy(t)|] < -at+e. 


veEV (A) 


Then the zero G-solution G[DES}, 
on subgraphs H; and Ho, but 


X y(t 
veV(HiU He) 


) 


Thus the zero G-solution G[DES}, 
subgraphs A; () H2. Similarly, let 


vEV (He) 


of differential equation system (DES},) is not sum-stable 


Ss 


S> Xp(t)]] + 


veEV(M1) 


Ss) X(t) 


vEV (A2) 


=€. 


of differential equation system (DES) is sum-stable on 
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I] xl < = and || S° X,(#)) <¢ 
veV(A1) veV (He) 

for a real number r. Then the zero G-solution G[DES},] of (DES1},) is not prod-stable on 

subgraphs H, and X,(t) is not bounded for v € V(Hp2) if r > 0. However, it is prod-stable on 


subgraphs AH; () H2 for 


I] *@l<) If xl! IT] x Ol =e 


veEV(AMiU 2) veV (Hi) veV (Hz) 


3.2 Linearized Differential Equations 


Applying these conclusions on linear differential equation systems in the previous section, we 
can find conditions on F;(X), 1 < i < m for the sum-stability and prod-stability at 0 following. 
For this objective, we need the following useful result. 


Lemma 3.5(([13]) Let X = AX + B(X) be a non-linear differential equation, where A is a 
constant n x n matric and Red; < 0 for all eigenvalues A; of A and B(X) is continuous defined 
ont >0, ||X|| < a with 

IB 


lim “</" — 9. 
|Ix\|0 || X| 


Then there exist constants c>0, 8 >0 and 6, 0<6<a such that 


|X (O)|| <e< implies that. ||X(t)|| < cee~®*/. 
Cc 


Theorem 3.6 Let (DES!) be a non-linear differential equation system, H a spanning subgraph 
of G[DES},] and 

F(X) = Fi (0) X+ R(X) 
such that 


|Ro(X)II 


lim us —9 
x0 |X| 


for tv € V(H). Then the zero G-solution of (DES},) is asymptotically sum-stable or asymp- 
totically prod-stable on H if Rea, <0 for each B,(t)e%* € By, v € V(H) in (DES). 


Proof Define c = max{c,, v € V(H)}, ¢ = min{e,, v © V(A)} and 6 = min{G,, v € 
V(H)}. Applying Lemma 3.5, we know that for Vu € V(A), 


é6 
\|X,(0)||<e< 2c implies that ||X,(t)|| < cee Bt/2, 
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Whence, 
S> X(t] < S> [X.()ll < |H]cce*? 
veV (A) veV(H) 
I] ©) < [TD 1Xe@]) < lel tler2, 
veV (A) veV(H) 
Consequently, 
lim S-> X,(t)|| +0 and lim I] x] -0. 
veV(H) veV(H) 


Thus the zero G-solution (DES?) is asymptotically sum-stable or asymptotically prod-stable 
on H by definition. 


3.3 Liapunov Functions on G-Solutions 


We have know Liapunov functions associated with differential equations. Similarly, we introduce 
Liapunov functions for determining the sum-stability or prod-stability of (DES},) following. 


Definition 3.7 Let (DES},) be a differential equation system, H < G[DES},] a spanning 
subgraph and a H-equilibrium point X* of (DES},). A differentiable function L: @ + R 
defined on an open subset OC C R” is called a Liapunov sum-function on X* for H if 


(1) L(X*) =0 ond =, *) >Oif YS X(t) # X*; 


vEV( veV(H) 
ail So x) <ofr CT KNHAX*, 
veV(H) veV(H) 
and a Liapunov prod-function on X* for H if 
(1) L(X*) =0 and L ( II 2) >Oif JT] X(t) 4x"; 
H) 


vEeV( veEV (HA) 


V(A) veV (A) 


a E( II (0) <0 fr I] Xv(t) 4 X*. 


Then, the following conclusions on the sum-stable and prod-stable of zero G-solutions of 
differential equations holds. 


Theorem 3.8 For a G-solution G[DES},| of a differential equation system (DES},) with 
initial value X,(0), let H be a spanning subgraph of G[DES},] and X* an equilibrium point of 
(DES}) on H. 


(1) If there is a Liapunov sum-function L : @ — R on X*, then the zero G-solution 
G[DES}] is sum-stable on X* for H. Furthermore, if 
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for So X(t) 4 X*, then the zero G-solution G[DES},] is asymptotically sum-stable on X* 
veEV (HA) 
for H. 


(2) If there is a Liapunov prod-function L: @ > R on X* for H, then the zero G-solution 
G[DES} |] is prod-stable on X* for H. Furthermore, if 


LY J] xX] <0 


veV (HA) 


for J] Xv(t) #4 X*, then the zero G-solution G[DES}] is asymptotically prod-stable on X* 
veEV (HA) 
for H. 


Proof Let « > 0 be aso small number that the closed ball B.(X*) centered at X* with 
radius ¢€ lies entirely in @ and w the minimum value of L on the boundary of B,(X“), i-e., 


the sphere S.(X*). Clearly, @ > 0 by assumption. Define U = {X € B,.(X*)|L(X) < a}. 
Notice that X* € U and L is non-increasing on > X,(t) by definition. Whence, there are 


veV(H) 
no solutions X,(t), v € V(H) starting in U such that S°> X,(t) meet the sphere S.(X%*). 
veEV(H) 
Thus all solutions X,(t), v € V(H) starting in U enable 5> X,(t) included in ball B,(X*). 
veEV(H) 


Consequently, the zero G-solution G[DES},] is sum-stable on H by definition. 


Now assume that 


El 3 a) | 20 


veV (HA) 
for So X,(t) # X*. Thus L is strictly decreasing on S> X,(t). If X,(t), ve V(A) are 
veEV(H) veV(H) 
solutions starting in U — X* such that > Xy(tn) — Y* for n — oo with Y* € B.(X%), 
veEV(H) 


then it must be Y* = X*. Otherwise, since 


by the assumption 
forall So X,(t) 4 X* and 
veV(H) 
Xy(tn) | > L(Y*) 


vEV 


by the continuity of L, if Y* # X*, let Y,(t),v © V(A) be the solutions starting at Y*. Then 


for any 7 > 0, 
LY > Mo) | SLO: 
veEV(H) 
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But then there is a contradiction 


L| So Xol(tn+n)) < LY") 
veEV (A) 
yields by letting Y,(0) =  S> Xuy(ty) for sufficiently large n. Thus, there must be Y,;* = X*. 
veV (HA) 
Whence, the zero G-solution G[DES},] is asymptotically sum-stable on H by definition. This 
is the conclusion (1). 


Similarly, we can prove the conclusion (2). 


The following result shows the combination of Liapunov sum-functions or prod-functions. 


Theorem 3.9 For a G-solution G[DES},| of a differential equation system (DES},) with 
initial value X_(0), let Hy, Hz be spanning subgraphs of G[DES},], X* an equilibrium point of 
(DES},) on Hi\J He and 

De a,j" y! 


i>0,j>0 
be a polynomial with a;,; > 0 for integers i,j > 0. Then Rt (Ly, L2) is a Liapunov sum-function 
or Liapunov prod-function on X* for Hy) H2 with conventions for integers i,j,k, > 0 that 


ay Li Ts S- X(t) | tawLPL} S- X,(t) 
veV (Hi U V(A2) veV (Hi U V(A2) 


=a;Li| > X,(t)| £4 
veV (A) oe 


veV (Hi) es 


if L1, Le are Liapunov sum-functions and 


ai; Li 13 II X,(t) | tally 
ve€V(A1U V(A2) 


=ajli| [J X.(t)| 24 
veEV (A) J. 


+anLh | S> X,(t) | £5 4 


+anlh} JT x @)25{ Jf x@ 


veV (Hi) veEV(H2) 


if L1,L2 are Liapunov prod-functions on X* for Hy, and Ho, respectively. Particularly, if 
there is a Liapunov sum-function (Liapunov prod-function) L on H, and Ho, then L is also a 


Liapunov sum-function (Liapunov prod-function) on Hy |) Ho. 
Proof Notice that 


dk = aij Geman: + jLATA Es) 
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if 74,7 > 1. Whence, 


ai; Li Lh S- X,(t)| >0 
veV (Ai U V (Ae) 
if 
Ii1{ S° X,(t)] 20 and Ig{ S> X(t)]} 0 
veV(A1) veV (2) 
and 
pegs 
Mo Tite) X,(t) ee 0 
veV (Ai U V (Ae) 

if 


Ei{| S> X(t)| <0 and bot S> X,(t)} <0. 
veEV (A) ve€V(H2) 
Thus R*(L 1, D2) is a Liapunov sum-function on X* for Hy LJ Ho. 
Similarly, we can know that Rt (Ly, L2) is a Liapunov prod-function on X* for Hy U He if 


[,, L2 are Liapunov prod-functions on X* for H, and Abo. 
Theorem 3.9 enables one easily to get the stability of the zero G-solutions of (DES}.). 


Corollary 3.10 For a differential equation system (DES), let H < G[DES}]] be a spanning 
subgraph. If Ly is a Liapunov function on vertex v for Vu € V(H), then the functions 


Lg= >> L, and LE = JT Ly 
veV(H) veV(H) 


are respectively Liapunov sum-function and Liapunov prod-function on graph H. Particularly, 
if L = Ly forVu € V(A), then L is both a Liapunov sum-function and a Liapunov prod-function 
on Hf. 


Example 3.11 Let (DES},) be determined by 


dx, /dt => Aqw21 da, /dt = A121 dx, /dt = Ani 21 
dix /dt = A1222 dx /dt = A22X2 dix / dt = An22 
diy, [dt = AinXn diy /dt = r2nXn diy /dt = XnnXn 


where all yj, 1 << i< m,1 <j < nare real and \y;, A Aij. if j1 A Jo for integers 1 <i <m. 
Let L= 27+ 03+---+22. Then 


L= Net + rin? ee rin 22, 


for integers 1 <i <n. Whence, it is a Liapunov function for the ith differential equation if 
rij < 0 for integers 1 < j < n. Now let H < G[LDES},] be a spanning subgraph of G[LDES}, ]. 
Then L is both a Liapunov sum-function and a Liapunov prod-function on H if A,; < 0 for 
Vu € V(#) by Corollaries 3.10. 
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Theorem 3.12 Let L: @ — R be a differentiable function with L(0) = 0 and L ( Sy x) > 
veEV(H) 


0 always holds in an area of its e-neighborhood U(e) of 0 for ¢ > 0, denoted by U*(0,e) such 


area of €-neighborhood of 0 with L SS X | >0andH < G[DES!] be a spanning subgraph. 
veEV(H) 


(1) If 
G S- X|||<™M 
veV(H) 


with M a positive number and 


i S> xX] >0 


veV(H) 


in U*(0,€), and for Ve > 0, there exists a positive number cy, ¢2 such that 


L . X |} >c, >0 implies L S- X |] >c>O0, 
veV(H) veV(H) 
then the zero G-solution G[DES},] is not sum-stable on H. Such a function L: @ > R is 


called a non-Liapunov sum-function on H. 


(2) If 


L| JJ x])<n 
veEV (HA) 
with N a positive number and 
L| [J x]>o 
vEV(H) 


in Ut (0,€), and for Ve > 0, there exists positive numbers d,, dz such that 


L II X | >d,>0 implies L II X | >d2>0, 
veV(H) veV(H) 
then the zero G-solution G[DES},] is not prod-stable on H. Such a function L: GC > R is 


called a non-Liapunov prod-function on H. 


Proof Generally, if ||L(X)|| is bounded and L (X) > 0 in U+(0,¢), and for Ve > 0, there 
exists positive numbers ci, c2 such that if D(X) > c > 0, then L (X) > cg > 0, we prove that 
there exists t; > to such that ||X(t1,to)|| > €o for a number €9 > 0, where X(t1, to) denotes 
the solution of (DES) passing through X (to). Otherwise, there must be ||X (£1, to)|| < €o for 
t > to. By L(X) > 0 we know that L(X(t)) > L(X(to)) > 0 for t > to. Combining this fact 
with the condition L(X) > cy > 0, we get that 

t 
LX (t)) = L(X(to)) + | FEE > r(e(t0)) + ealt = to) 


to 
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Thus L(X (t)) too if t Foo, a contradiction to the assumption that L(X) is bounded. 
Whence, there exists t; > to such that 


|| X(t1, to) || > €0. 


Applying this conclusion, we immediately know that the zero G-solution G[DES}]] is not sum- 


stable or prod-stable on H by conditions in (1) or (2). 


Similar to Theorem 3.9, we know results for non-Liapunov sum-function or prod-function 
by Theorem 3.12 following. 


Theorem 3.13 For a G-solution G[DES},] of a differential equation system (DES},) with 
initial value X,(0), let Hy, Hz be spanning subgraphs of G[DES},], 0 an equilibrium point of 
(DES!) on Hi He. Then Rt(L1,L2) is a non-Liapunov sum-function or non-Liapunov 
prod-function on 0 for Hi \) Hz with conventions for 


ajl te S- RAD ag, S- X,(t) 
v€V (M1 U V(A2) veV (MU V(A2) 
and 
ai; Li 3 II RG) ag ler, II X,(t) 
v€V (M1 U V(A2) veV (MU V(A2) 


the same as in Theorem 3.9 if L1, Lo are non-Liapunov sum-functions or non-Liapunov prod- 
functions on 0 for Hy and Hz, respectively. Particularly, if there is a non-Liapunov sum- 
function (non-Liapunov prod-function) L on Hy and Ho, then L is also a non-Liapunov sum- 


function (non-Liapunov prod-function) on Hy VU He. 


Proof Similarly, we can show that R*(L , D2) satisfies these conditions on HU H2 for 
non-Liapunov sum-functions or non-Liapunov prod-functions in Theorem 3.12 if [1,2 are 
non-Liapunov sum-functions or non-Liapunov prod-functions on 0 for H; and He, respectively. 


Thus Rt (Ly, L2) is a non-Liapunov sum-function or non-Liapunov prod-function on 0. 


Corollary 3.14 For a differential equation system (DES},), let H < G[DES}]] be a spanning 
subgraph. If Ly is a non-Liapunov function on vertex v for Vu € V(H), then the functions 


LE= S° L, and LB = JI L, 


veV(H) veV(H) 


are respectively non-Liapunov sum-function and non-Liapunov prod-function on graph H. Par- 
ticularly, if L = Ly for Vv € V(X), then L is both a non-Liapunov sum-function and a non- 


Liapunov prod-function on H. 


Example 3.15 Let (DES!) be 


. 2 2 . 2 2 . 2 2 
ty = A\yay — A125 te = Aguz — A2x5 t1 = Am@{— Am®5 
r : x eee ; r 


a uy 2 m 
t2 = Zt12%2 2 = Zt12%2 LQ = VW F142 


2 2 2 
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with constants \; > 0 for integers 1 < i < m and L(a1,22) = 2? — 242. Then L(21, 22) = 
4d;2,L (21,22) for the i-th equation in (DES},). Calculation shows that L(a1,x72) > 0 if 
ay > V2x2 or 41 < —V2x2 and L(21, £2) > 4c? for L(x1,2%2) > cin the area of L(x, x2) > 0. 
Applying Theorem 3.12, we know the zero solution of (DES},) is not stable for the i-th equation 
for any integer 1 <i < m. Applying Corollary 3.14, we know that L is a non-Liapunov sum- 
function and non-Liapunov prod-function on any spanning subgraph H < G[DES}]]. 


§4. Global Stability of Shifted Non-Solvable Differential Equations 


The differential equation systems (DES},) discussed in previous sections are all in a same 
Euclidean space R”. We consider the case that they are not in a same space R”, i.e., shifted 
differential equation systems in this section. These differential equation systems and their 
non-solvability are defined in the following. 


Definition 4.1 A shifted differential equation system (SDES},) is such a differential equation 


system 


X, = F,(X1), Xo = Fo(X2),-++ Xm = Fn(Xm) (SDES}) 
with 

Xy = (21, %2,°°* , 81, %10141), £1042))*°* Lin), 

X2q = (&1,£2,°°* , £1, F2(141), La142),°** » Lan), 

Xm = (£1, %2,°°* » L1, Xm (141)> Pm(l4+2)9°°* :tmn); 
where £1, %2,°*+ , 1, X45), 1 <i <m,1<j <n-—l are distinct variables and F, :R” > R” 


is continuous such that F,(0) =0 for integers 1<s<m. 

A shifted differential equation system (SDES},) is non-solvable if there are integers i,j, 1 < 
i,j <m and an integer k, 1 <k <1 such that oll (t) # oll] (t), where oll (t), wll] (t) are solutions 
xp(t) of the i-th and j-th equations in (SDES},), respectively. 


The number dim(SDESj,) of variables #1, %2,-++ ,@1,Tii4j),1 <i<mil<j<n-lin 
Definition 4.1 is uniquely determined by (SDES1), ie., dim(SDES},) = mn —(m-— 1)l. For 
classifying and finding the stability of these differential equations, we similarly introduce the 
linearized basis graphs G[S DES] of a shifted differential equation system to that of (DES}), 
i.e., a vertex-edge labeled graph with 


V(G|SDES1]) ={B|l<i<m}, 
B(G[SDES}]) = {(Bi, B)| BOB; #V1<S 4,5 < mj, 


where &; is the solution basis of the i-th linearized differential equation X; = F/(0)X; for 
integers 1 < i < m, called such a vertex-edge labeled graph G[SDES},] the G-solution of 
(SDES},) and its zero G-solution replaced 4; by (0,--- ,0) (|B;| times) and (1) AB; by 
(0,--- ,0) (|Ai(] B;| times) for integers 1 < i,j < m. 

Let (LDES}), (LDES} )' be linearized differential equation systems of shifted differential 
equation systems (SDES}) and (SDES},) with G-solutions H, H’. Similarly, they are called 
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combinatorially equivalent if there is an isomorphism y : H — H’ of graph and labelings 
6, 7 on H and H’ respectively such that y6(~) = Ty(x) for Vz € V(H)U E(A), denoted by 
(SDES},) © (SDES1,)'. Notice that if we remove these superfluous variables from G[SDES}], 
then we get nothing but the same vertex-edge labeled graph of (LDES},) in R!. Thus we can 
classify shifted differential similarly to (LDES},) in R!. The following result can be proved 
similarly to Theorem 2.14. 


Theorem 4.2. Let (LDES},), (LDES,,)' be linearized differential equation systems of two 
shifted differential equation systems (SDES},), (SDES1,)' with integral labeled graphs H, H’. 
Then (SDES1) & (SDES1,) if and only if H = H’. 


The stability of these shifted differential equation systems (S DES} ) is also similarly to 
that of (DES},). For example, we know the results on the stability of (SDES},) similar to 
Theorems 2.22, 2.27 and 3.6 following. 


Theorem 4.3 Let (LDES},) be a shifted linear differential equation systems and H < G[LDES}] 
a spanning subgraph. A zero G-solution of (LDES},) is asymptotically sum-stable on H if and 
only if Rea, <0 for each 3,,(t)e%* € B, in (LDES") hold for Vv € V(H) and it is asymptot- 
ically prod-stable on H if and only if > Reay <0 for each B,(t)e%! € B, in (LDES?). 
veV(H) 
Theorem 4.4 Let (SDES},) be a shifted differential equation system, H < G[SDES},] a 
spanning subgraph and 
F(X) = F, (0) X + R,(X) 
such that Ry(X)| 7 
WxI—o ||| 
for Vv € V(H). Then the zero G-solution of (SDES},) is asymptotically sum-stable or asymp- 
totically prod-stable on H if Rea, <0 for each 3,,(t)e%* € By, v€ V(H) in (SDES}). 


For the Liapunov sum-function or Liapunov prod-function of a shifted differential equation 
system (SDES}_), we choose it to be a differentiable function L : @ C Réim(SDES;,) _, R with 
conditions in Definition 3.7 hold. Then we know the following result similar to Theorem 3.8. 


Theorem 4.5 For a G-solution G[SDES})] of a shifted differential equation system (SDES} ) 
with initial value X,(0), let H be a spanning subgraph of G[DES},] and X* an equilibrium 
point of (SDES}) on H. 


(1) If there is a Liapunov sum-function L: @ Cc Riim(SDES;,) _, Ron X*, then the zero 
G-solution G[SDES}] is sum-stable on X* for H, and furthermore, if 


7) i eee can ka 


veV (HA) 


for SY > X(t) # X*, then the zero G-solution G[SDES},| is asymptotically sum-stable on 
veV(H) 
X* for H. 
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(2) If there is a Liapunov prod-function L : @ Cc Rim(SDES,,)  R on X* for H, then 
the zero G-solution G[SDES}] is prod-stable on X* for H, and furthermore, if 


LY JJ X.()} <0 
veEV (HA) 
for |] Xv(t) 4 X*, then the zero G-solution G[SDES}] is asymptotically prod-stable on 


veEV (HA) 
X* for H. 


§5. Applications 


5.1 Global Control of Infectious Diseases 


An immediate application of non-solvable differential equations is the globally control of infec- 
tious diseases with more than one infectious virus in an area. Assume that there are three kind 
groups in persons at time ¢, ie., infected I(t), susceptible S(t) and recovered R(t), and the 
total population is constant in that area. We consider two cases of virus for infectious diseases: 


Case 1 There are m known virus \,, 42,°-+,Vm with infected rate k;, heal rate h; for integers 


1<i<m and an person infected a virus V; will never infects other viruses V; for 7 #4. 


Case 2 There are m varying V\,,%2,--: ,Vm from a virus V with infected rate k;, heal rate hy 


for integers 1<%1<m such as those shown in Fig.5.1. 


Q-®---@ 


Fig.5.1 


We are easily to establish a non-solvable differential model for the spread of infectious 
viruses by applying the SIR model of one infectious disease following: 


S$ =-k,SI S$ =—k)SI S = —kpSI 
I =k SI —hil T=kSI-hoI °-:: I = kmSI — mI (DES}) 
R=hl R= hol R=h,I 


Notice that the total population is constant by assumption, i.e., S++ R is constant. 
Thus we only need to consider the following simplified system 


S=-k SI S=—-kSI S =—kmSI 


; , (DES) 
f=hSI-hyI T = kySI — hel T= kp SI — hmI 


The equilibrium points of this system are I = 0, the S-axis with linearization at equilibrium 
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points 


S=-kS S=—-hkS S =—kmS 
; ve (LDES),) 
T=khS—h, T=kS — he T=knS—hn 


Calculation shows that the eigenvalues of the 7th equation are 0 and k;S—h;, which is negative, 
i.e., stable if 0 < S < h;/k; for integers 1 <i < m. For any spanning subgraph H < G[LDES}], 
we know that its zero G-solution is asymptotically sum-stable on H if 0 < S < hy/ky for 
v € V(H) by Theorem 2.22, and it is asymptotically sum-stable on H if 


Y) (S-m)<0 ie, O<s< Yom / YO 


veEV (HA) veEV (A) veV(H) 
by Theorem 2.27. Notice that if [;(t), S;(t) are probability functions for infectious viruses 
¥,, 1 <i < min an area, then [] Z(t) and [] 5S;(t) are just the probability functions for 


i=1 = 
all these infectious viruses. This fact enables one to get the conclusion following for globally 


control of infectious diseases. 


Conclusion 5.1 For m infectious viruses V,,%2,--- ,¥%m in an area with infected rate k;, heal 
rate hy for integers 1 <i<m, then they decline to 0 finally if 


0<8<S0h Se 
1=1 t=1 


i.e., these infectious viruses are globally controlled. Particularly, they are globally controlled if 
each of them is controlled in this area. 


5.2 Dynamical Equations of Instable Structure 


There are two kind of engineering structures, i.e., stable and instable. An engineering structure 
is instable if its state moving further away and the equilibrium is upset after being moved 
slightly. For example, the structure (a) is engineering stable but (b) is not shown in Fig.5.2, 


A @ 


Fig.5.2 


where each edge is a rigid body and each vertex denotes a hinged connection. The motion of 
a stable structure can be characterized similarly as a rigid body. But such a way can not be 


applied for instable structures for their internal deformations such as those shown in Fig.5.3. 
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Fig.5.3 


Furthermore, let Y,, Yo,--: ,Am be m particles in R° with some relations, for instance, 
the gravitation between particles Y; and #; for 1 < i, 7 < m. Thus we get an instable structure 
underlying a graph G with 

V(G) ={F1, Pr,-++, Pm}; 

E(G) = {(A%, Y;)|there exists a relation between Y; and Y;}. 

For example, the underlying graph in Fig.5.4 is Cy. Assume the dynamical behavior of particle 


P, at time t has been completely characterized by the differential equations X = F,(X,t), 
where X = (21,%2,23). Then we get a non-solvable differential equation system 


underlying the graph G. Particularly, if all differential equations are autonomous, i.e., depend 
on X alone, not on time t, we get a non-solvable autonomous differential equation system 


X=F,(X), 1<i<m. 


All of these differential equation systems particularly answer a question presented in [3] for 
establishing the graph dynamics, and if they satisfy conditions in Theorems 2.22, 2.27 or 3.6, 
then they are sum-stable or prod-stable. For example, let the motion equations of 4 members 
in Fig.5.3 be respectively 


AB: Xap = 0; CD: Xcp = 0, AC: Xac = aac, BC: Xgc =agc, 


where X4p,Xcp,Xac and Xgc denote central positions of members AB,CD, AC, BC and 
aAc,@pc are constants. Solving these equations enable one to get 


Xap =capt+dap, Xac =aact? +cact+ dac, 


Xcp =ccpt+dep, Xsc =azct? + cect + dzc, 


where c4p,cac,ccp,cBc,¢d4B,d4c,dcp,dgc are constants. Thus we get a non-solvable dif- 
ferential equation system 


X=0; X=0, X=aac, X =azc, 
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or a non-solvable algebraic equation system 


X =capt+dap, X =a,ct? + cact+dac, 


X=copt+dcp, X= apct? +cpct+ dsc 


for characterizing the behavior of the instable structure in Fig.5.3 if constants c4p,CAc, CCD; CBC; 
dap,dac,dcp,dgc are different. 
Now let X1, X2,--- , Xm be the respectively positions in R? with initial values X?, X9,--- , X° 


X9, X9,---, X° and Mj, Mo,--- , Mm the masses of particles P:, Po,--- , Pm. Ifm = 2, then 


from Newton’s law of gravitation we get that 


where G is the gravitational constant. Let X = X2— X1 = (#1, 22,23). Calculation shows that 


é Xx 
X = —G(M,+4+ M2) 3° 
|X| 


Such an equation can be completely solved by introducing the spherical polar coordinates 


1 =rcos@cosé 
r= rcos@cosé 


23 =rsind 


with r > 0,0 < 6 < 7,0 < 0 < 27, where r = ||X||, 6 = 2X0z,0 = ZX'ox with X' 
the projection of X in the plane xoy are parameters with r = a/(1 + ecos@) hold for some 


constants a,¢. Whence, 


x(0) Gah [ ( rae it and X(t) =-can, f ( | rae a 


Notice the additivity of gravitation between particles. The gravitational action of particles 
Py, Po,+++, Am on FY can be regarded as the respective actions of Py, Wo,--- , Am on FY, 
such as those shown in Fig.5.4. 


Fig.5.4 
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Thus we can establish the differential equations two by two, i.e., A, acts on Y, Wg acts on 


P,-++, Am acts on Y and get a non-solvable differential equation system 
: Sea a 
X =GM;——_,,  B#PY, 1<i<m. 
[Xi — X| 


Fortunately, each of these differential equations in this system can be solved likewise that of 
m = 2. Not loss of generality, assume X(t) to be the solution of the differential equation in 
the case of A; 4 YP, 1<i<m. Then 


a XxX; —X 
xo= % R= Yu f (/ Aen) 4 
PAP PAP |Xi — X| 
is nothing but the position of particle Y at time t in R? under the actions of A, 4 Y for 
integers 1 < i < m, i.e., its position can be characterized completely by the additivity of 


gravitational force. 
5.3 Global Stability of Multilateral Matters 


Usually, one determines the behavior of a matter by observing its appearances revealed before 
one’s eyes. If a matter emerges more lateralities before one’s eyes, for instance the different 
states of a multiple state matter. We have to establish different models, particularly, differential 
equations for understanding that matter. In fact, each of these differential equations can be 
solved but they are contradictory altogether, i-e., non-solvable in common meaning. Such a 
multilateral matter is globally stable if these differential equations are sum or prod-stable in all. 

Concretely, let S$ ,52,---, 5m be m lateral appearances of a matter .W in R® which are 
respectively characterized by differential equations 


where X; € R°, a 3-dimensional vector of surveying parameters for S;,1 <i< m. Thus we get 


a non-solvable differential equations 


in R3. Noticing that all these equations characterize a same matter .@, there must be equilib- 
rium points X* for all these equations. Let 


H,(X,t) = Hi(X*)X + R,(X*), 


where 
ny al ny 
A 


is an n X n matrix. Consider the non-solvable linear differential equation system 


X =Hi(X*)X, 1<i<m (LDES! ) 
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with a basis graph G. According to Theorem 3.6, if 


— IRQOI 
ix xn 


for integers 1 < i < m, then the G-solution of these differential equations is asymptotically 
sum-stable or asymptotically prod-stable on G if each Real! < 0 for all eigenvalues all of 
matrix H/(X*), 1<i<m. Thus we therefore determine the behavior of matter @ is globally 
stable nearly enough X*. Otherwise, if there exists such an equation which is not stable at the 
point X*, then the matter WH is not globally stable. By such a way, if we can determine these 
differential equations are stable in everywhere, then we can finally conclude that M is globally 
stable. 

Conversely, let @ be a globally stable matter characterized by a non-solvable differential 
equation 

X = H,(X,t) 


for its laterality S;, 1<7<m. Then the differential equations 
X=H,(X,t), 1l<i<m (DES}) 


are sum-stable or prod-stable in all by definition. Consequently, we get a sum-stable or prod- 
stable non-solvable differential equation system. 


Combining all of these previous discussions, we get an interesting conclusion following. 


Conclusion 5.2 Let ME A” be respectively the sets of globally stable multilateral matters, 
non-stable multilateral matters characterized by non-solvable differential equation systems and 
QGE,DVE the sets of sum or prod-stable non-solvable differential equation systems, not sum or 
prod-stable non-solvable differential equation systems. then 

(1) V4 6 AE = A(DES}) € GE; 

(2) VM eM” > DES!) € GE. 

Particularly, let @ be a multiple state matter. If all of its states are stable, then @ is 
globally stable. Otherwise, it is unstable. 
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